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One of the most direct way to constrain the matter density of the universe 
Dm is to measure the baryon content in X-Rays clusters of galaxies. Typical 
value of the mean gas mass fraction is / gas ~ 0.15/V 3 / 2 which leads to O m < 

0.4. In this talk I will discuss the issue of the gas fraction radial distribution 
inside a cluster and foccus on the apparent discrepancy between the theoretical 
shape predicted by numerical simulations and the observations. I will show 
that a significant part of this discrepancy is due to several systematics in both 
the gas mass and total mass determinations. Revising the gas as well as the 
binding masses removes such a discrepancy and results in a lower baryon mass 
fraction of / b (r 500 ) ~ 10% (Sadat & Blanchard 2001). Finally, I will discuss 
the cosmological constraints we obtain when we combine the revised baryon 
fraction with primordial fluctuations. 
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1 Introduction 


Clusters of galaxies are the most largest virialized systems in the Universe, hence the 
measurement of the baryonic to total mass ratio is representative of the cosmic value 

Qbbn 

fb = - 77 —. The combination of the universal baryon density as predicted from light 
element abundances through the theory of big-bang nucleosynthesis with the cluster gas 
mass fraction, lead to a direct determination of f2 m . This test to is in principle the 
most direct one as it depends on very few reasonable assumptions: the cluster formation 
by collapse from a well mixed medium, and no segregation between the gas and the dark 
matter as shown by numerical simulations. 

In their pioneering paper, White and his collaborators have used Rosat observations of 
Coma cluster to estimate the gas mass fraction they found that the archetypical coma 
cluster has a baryonic content of Mf, ~ 3 . 4 h 5 0 5 ^ 2 10 14 M Q with a very negligible contribution 
from stars iff* ~ 2 . 2 /i 5 0 1 10 13 Mq and a binding mass of M to t ~ 2 . 2 /i 5 0 1 10 15 Mq. Hence they 
derived a gas fraction of / gas ~ 0.15 h§^ 2 at Abell radius, which is several times higher 
than expected from the observed light element abundances (0.024 < fib ^Jo A 0.08) in 
Q rn = 1 universe, leading to the conclusion that we live in a low density universe (White et 
al. 1993). Several analyses of large samples of clusters have confirmed such a high baryon 
fraction with values ranging from 10 — 20%. These determinations depend on the radius 
where the gas mass fraction is estimated, since the gas fraction slightly increases with 
radius, but even when the gas fraction is determined in the same equivalent radius say r 50 o 
(radius encompassing a mean density 500 times the critical density), current estimations 
are still discrepant (Evrard 1997). It has also be argued that the baryon fraction increases 
with cluster temperature (mass) (David et al. 1995, Mohr et al. 1999, Arnaud & Evrard 
1999). This has often been explained as due to non-gravitational processes such as heating 
effects occuring during galaxy cluster formation. In this talk I will address the issue of 
the baryon fraction in clusters by examining critically its distribution within the cluster 
and show how some systematic effects can bias the gas mass fraction determination. In 
the first section I will review the classical analysis method used to estimate gas mass 
fractions. Comparison between numerical simulations and observation is provided in 
section 2. Section 3 presents the systematic uncertainties on various steps involved in 
the gas mass fraction derivation and section 4 describes the revised gas mass fraction 
and the consequence on both the shape of the gas mass fraction profile as well as on 
the mean mass density of the universe. In section 5 we investigate the constraints on 
cosmological parameters by combining the cosmic microwave background (CMB) with 
the galaxy cluster baryon fraction. 
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2 Gas mass fraction determination 


The standard way to derive the gas fraction proceeds from the following steps: 

• Assume a (3 — model to fit the observed surface brightness Sj^l + (#/# c ) 2 ) -3/3+1 / 2 where 
6 is the projected angular distance to the center. The gas density profile (hence the gas 
mass profile) is then derived using the two best-fit parameters (3 and 9 C . 

• Use the hydrostatic equation to derive the total mass inside the same radius: M tot (r) = 

elm. PTyj' + ' The t°t a I mass th us depends linearly on both (3 and T x . 

Hence, if the slope of the gas density is poorly determined, it will have a drastic influence 
on the derived mass. 


2.1 Shortcomings of this method 

The X-ray gas emission is rarely traced out to very large radii, so that gas fractions 
are measured only up to an X-rays limiting radius for which the signal-to-noise is good 
enough. Most of the time, it was necessary to extrapolate to the virial radius r 2 oo or to 
some other outer radius like r 50 o- It is then essential to realize that gas mass is estimated 
at a radius at which the actual emission is poorly constrained, with a value of the order 
of the X-rays background or less. Moreover, the emissivity being dominated by the cen¬ 
tral region, a parametric fit will be rather insensitive to the outer part of the gas profile. 
Finally, gas masses are estimated assuming that the gas is smoothly distributed within 
the intracluster medium. Concerning total mass determination, the /3 - model approach 
introduces an intrinsic scatter in the total mass estimates, of ~ 30% at r 500 which can be 
considerably reduced by using scaling relations between the mass and the temperature T x 
calibrated with numerical simulations (Evrard 1998). Moreover, the f3 - model approach 
leads to total masses which are systematically lower than what one would obtained using 
scaling relations between the mass and the temperature T x calibrated with numerical 
simulations (Roussel et al. 2000, hereafter RSB00). However, some problems exist in 
this latter approach. There is a scatter in the M — T x relation existing among different 
authors: T x = 4.75(M 200 /10 15 M Q ) 2 / 3 keV (Evrard, Metzler & Navarro 1996, hereafter 
EMN96) and T x = 3.81(Mi7 8 /lO 15 M 0 ) 2 / 3 keV (Bryan and Norman 1998, BN98); these 
two normalizations being the most extreme values. There is also the problem that most 
of the numerical simulations do not include the gas cooling neither the heating of the gas 
with supernovae or AGN, which may introduce a further change in the relation, although 
these effects might play minor role in the case of hot and massive clusters. 
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3 The baryon fraction profile: Observations versus 
simulations 

In the case where only gravity is acting during cluster formation, the gas fraction profile 
is expected to follow a scaling law, depending only on the contrast density 5. Hydro- 
dynamical simulations of an X-rays cluster from different groups (the Santa Barbara 
group) have shown that the gas fraction normalized to the global value, increases in the 
inner part and then tends to flatten in the outer part to reach the cosmic value see Fig. 
[T]. On the observational side, the distribution of the gas fraction within clusters has been 
widely examined by RSBOO for a large sample of groups and galaxy clusters. The result¬ 
ing shape is shown in (Fig. |lj). To ensure consistency gas fractions - calibrated by the 
numerical simulation - published in RSBOO are used. For comparison gas fractions at r 50 o 
and r 2 oo from other authors are also plotted in (Fig. |IJ). The comparison of both profiles 
is very surprising as it shows that: the gas fraction profile derived from observations is in 
strong disagreement with the numerical simulations results with a global value of 16%, its 
shape continuously increases and does not exhibit the flattening in the outer parts as seen 
in numerical simulations. Clearly, this discrepancy calls for caution when one is using the 
gas fraction to set upper limit on the mean density of the universe. 

One possibility is that processes acting during the cluster formation are not well under¬ 
stood and probably non-gravitational processes could have played an important role. But 
this is probably not the case, because large departure from scaling laws is expected which 
is not observed in the data (RSBOO). 


4 Uncertainties in gas mass estimation 

4.1 The extrapolation problem 

X-rays gas masses are usually derived up to an X-rays limiting radius Rxiim where the 
signal-to-noise is good enough. Most of the time it was necessay to extrapolate the 
/3-model up to the virial radius r 2 oo, where the emission is actually dominated by the 
background. Moreover a parametric fit will be rather insensitive to the outer part, because 
the emissivity is dominated by the central region. This crude extrapolation may then raise 
some question on its validity. In a recent study of a large sample of X-rays ROSAT images 
in which the data trace the emission up to very large radii (approximately up to the virial 
radius), it has been shown that the /3-model does not provide an accurate description 
of the surface brightness over the whole range of radii. They have found that the outer 
slopes (0.3r 20 o to r 20 o) are actually steeper than the slope found when the /3-model is 
applied to the whole cluster. Hence, their derived gas masses are lower than what one 
would obtain by extrapolating the /3-model. 


3 


0.20 



Figure 1: The distribution of the gas fraction f g versus 6 (5 = p(< r)/p c where p c = 
3H 2 /8nG). The observations are from RSBOO filled circles, Ettori&Fabian 1999 empty 
triangle and Arnaud & Evrard 1999 crosses. Statistical uncertainties on these quantities 
are smaller than the symbols size. The squares connected by a dashed line correspond to 
the theoretical gas fraction for a global value of 16%. Small squares connected by a line 
correspond to the gas fraction calculated from numerical simulations including winds by 
Metzler & Evrard (1997): for a 6 keV cluster ( long dashes ), a 3 keV cluster ( small dashes ) 
and without including winds ( continuous line) with a global value of 20 %. 

4.2 The clumping effect 

In most of the studies, the ICM gas is assumed to be uniform, which is not true, as actual 
clusters do exhibit a certain level of density fluctuations at small and large scale due 
probably to accretion and mergers events. This clumping is assumed to be more pronouced 
in the outer parts of the clusters where the relaxation is not completely achieved. In the 
presence of clumping the gas masses are overestimated by a factor \[C where C =< 
p 2 g > / < p g > 2 (p g is the gas density). This factor has been quantified from numerical 
simulations by Mathiesen et al. (1999). They found that models assuming a uniform 
density introduce a bias of a [C ~ 1.16 at S = 500. 
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5 The mean gas fraction revisited: consequence on 
gas fraction profile and 

Here, the gas masses are corrected for the extrapolation problem by using the Vikhlinin 
et ah (1999) sample which is the only sample for which the emissivity is detected up 
to virial regions and then corrected the gas mass at S = 500 for the clumping effect 
using the Mathiesen et al. (1999) factor. The total mass is derived from scaling relations 
calibrated from numerical simulations (using both EMN96 and BN98 calibrations). The 
result is plotted in figure |2|. From this figure, we can see that the gas fraction distribution, 
in particular the asymptotic behaviour, is now consistent with numerical simulations 
predictions, (although a slight difference remains in the inner parts at 6 < 10 4 which 
can well be due to energy injection, not negligible in these regions). The resulting gas 
fraction at r 50 o is now lower than previous estimates. By matching the data point at 
5 = 500, we find / gas = 0.0875 ± 0.0075 (0.108 ± 0.0092) at 68% confidence level with 
BN98 (EMN96) normalization, the uncertainty being due to numerical simulations (Fig. 
|2|). Adding the stellar mean contribution of 1%, our final value of the baryon fraction at 
the virial radius is: % ~ 0.10 ± 0.01 (with BN98 calibration) and % ~ 0.12 ± 0.01 (with 
EMN96 calibration). Combining the derived baryon fraction with the universal baryon 
density 1% ~ 0.08 as derived from recent D/H measurement reported by O’Meara et al. 
2000 (this value has been confirmed by recent results on CMB fluctuations measurements 
from Boomerang mission 2001), we find the following value for fio : Ho = 0.8 ± 0.1 using 
the baryon fraction estimated with BN98 normalization. 


6 Combining CMB with the baryon fraction: new 
constraints on the cosmological parameters 

Given the fact that CMB implies degeneracies between the fundamental cosmological 
parameters Hq, Ha and fi m , (see Fig. 3) it is interesting to look for combination with 
constraint that does not measure directly any of these parameters. Here 1 will present the 
results we obtain when we investigate the combination with the baryon fraction inferred 
from cluster observations using both the high baryon fraction as found in the literature 
and the revised low baryon fraction inferred in this work (Douspis et al. 2001). 

6.0.1 The high baryon fraction case 

The main result from adding the constraint from the baryon fraction is to restrict signif¬ 
icantly the various contours obtained with CMB alone (see figure ??). This is very clear 
from the Ho — r /io plane: the area of the contours is much reduced, leading to a preferred 
model which is around Ho ~ 60 km/s/Mpc and 7710 ~ 5. The other two-dimensional plots 
also reveal that all the contours are significantly reduced: actually most of the one sigma 
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Figure 2: The distribution of the gas fraction f g versus 5: circles correspond to the 
observed f g from RSBOO but restricted to the region where the gas is detected (no extrap¬ 
olation), using EMN96 (empty circles ) and BN98 (filled circles ) calibration. Empty ( filled ) 
star corresponds to the estimated f g using VFJ99 sample and corrected for the clumping 
with EMN96 (BN98) calibration. Filled squares linked by a dashed line correspond to the 
gas fraction predicted by numerical simulations for a global fraction of 8.75%. The upper 
dashed line corresponds to the case where the global fraction is 10 . 8 % 
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contours are of the order of the grid size so that the best model and la ranges should 
be used with caution. It is highly interesting that the cosmological model is now very 
well constrained: fR ~ 0.6, Q to t ~ 1- (implying VL m ~ 0.4, H 0 ~ 60 km/s/Mpc, 7710 ~ 5, 
n ~ 0.85). As we mentioned the allowed range have to be interpreted with caution, yet it 
is worth noticing that the preferred value of the Hubble constant is rather low and that 
higher values like 80 km/s/Mpc lie uncomfortably outside of the preferred region. 

6.0.2 The low baryon fraction case 

Using the low value for the baryon fraction leads similarly to much restricted contours. 
However the latter differ from the contours obtained with a higher baryon fraction : the 
preferred value for r] w is again close to 5., the preferred value of the Hubble constant is now 
very low (H 0 ~ 40 km/s/Mpc). It is worth noticing that values as H 0 ~ 60 km/s/Mpc 
are still at the edge of the 3 sigma contour and for this reason cannot be entirely ruled 
out. The preferred model is ~ 0-3, 1-hot ~ 1.1. This implies an Q m ~ 0.8 with H 0 ~ 40 
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Figure 3: Different contour plots from the analysis of the most recent CMB data (COBE, 
BOOMERANG, MAXIMA, DASI). The dashed red line define the 68, 95 and 99% GCL 
when projected on the axis. The blue filled contours are the corresponding confidence 
intervals in 2 dimensions. 


km/s/Mpc, 7/10 ~ 5. and n ~ 0.85 (see figure 4). 


7 Conclusions 

The baryon fraction in X-rays clusters is one of the most direct way to constrain the 
density of the universe. Previous estimations of the baryon fraction yield to a high value 
/b ~ 15% h 50 3 ^ 2 , favouring low density universe Do ~ 0.3. In my talk I have shown that 
the radial prohl of the gas fraction as derived from observations differs strongly from 
numerical simulations results. I argue this is not due to non-gravitational processes which 
may occur during the cluster formation but is rather due to various systematics such as 
the extrapolation of a /3-model and the dumpiness of the gas that have not been taken into 
account in previous estimations. Correcting the gas masses from these biases and using 
total masses calibrated from numerical simulations yield to a better agreement between 
observations and numerical simulations for a global value of the order of rs j io(i2)% fc 5 y /2 
depending on which normalization we use to compute the total mass. This revised value 
is lower than previous ones and leads to Do r^j 0.8(0.66) h 5 - 0 1/2 . 

Finally combining the most recent CMB data with the baryon fraction leads to very 


7 














Ho 




Ho Ho 



Ho 






A Ho A H 0 


Figure 4: Combined analysis of CMB and the baryon fraction. Left: the high baryon 
fraction value case with /b « 15%, right : the low baryon fraction case with /b ~ 10%. 
The definition of contours is same as in Fig. 3 

interesting constraints on the various parameters. Using the revised value (the low baryon 
fraction value), we found that the best model is fR ~ 0.3, tt tot ~ 1.1 , H 0 ~ 40 km/s/Mpc, 
rji o ~ 5., n ~ 0.85 implying a high density parameter fl m ~ 0.8 consistent with the 
determination from the baryon fraction (Sadat & Blanchard 2001) and cluster abundance 
evolution (Blanchard et al. 2000). 
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